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I.  INTRODUCTION 


We  have  been  concerned  with  the  modeling  of  premixed,  laminar,  steady- 
state  flames.  1  In  a  previous  paper, ^  hereinafter  called  paper  I,  we  studied  a 
set  of  five  H2-02-N2  flames  and  found  that  the  species  profiles  were  very 
nearly  independent  of  the  transport  algorithms  used  to  describe  the 
multispecies  transport  phenomena.  The  computed  flame  speeds  corresponding  to 
this  set  span  a  range  of  more  than  a  factor  of  60  in  magnitude.  Because  no 
exceptions  could  be  found  we  surmised  that  this  independence  of  transport 
algorithm  (not  of  transport  parameters)  was  of  a  general  nature.  We  wished 
to  verify  this  inference  by  performing  a  similar  computer  study  for  the 
methane-air  flame.  The  results  are  reported  here. 

Methane-air  was  selected  because  there  is  sufficient  information  to  model 
it  and  because  the  species  therein  form  a  significant  subset  of  the  modeling 
of  higher  hydrocarbon  flames.*^  The  method  of  approach  is  to  select  and  fix 
input  parameters  for  CH^-air  mixtures  and  then  vary  the  mathematical 
approximations  to  the  multicomponent  transport  properties.  We  consider  only 
flames  whose  total  pressure  is  one  atmosphere  and  whose  unburned  gas 
temperature  is  298K. 

In  paper  I  we  employed  six  methods  of  approximating  transport.  Five  of 
these  methods  encompassed  the  range  of  approximations  commonly  found  in  the 
literature.  The  sixth  is  a  method  we  developed  that  is  computationally 
efficient  and  provides  results  extremely  close  to  the  most  precise  method. 
Because  of  our  previous  experience,  (paper  I)  we  elected  here  to  use  only  the 
most  precise,  the  most  simple  and  our  best  compromise  transport  methods.  They 
are,  in  the  notation  of  paper  I:  Method  I,  Method  V  and  Method  VI, 
respectively. 

To  provide  as  great  a  range  of  flame  speeds  as  possible,  we  have  also 
examined  the  stoichiometric  methane-oxygen  system  in  buffer  gases  of  argon  and 
helium. 

The  remainder  of  this  paper  consists  of  four  sections.  Section  II 
briefly  recalls  the  appropriate  flame  equations.  Section  III  lists  the  input 
parameters  used.  Section  IV  reviews  the  approximate  transport  methods,  and 
the  results  and  discussion  are  found  in  Section  V. 


^J.M.  Heimevl  and  T.P.  Coffee,  " The  Detailed  Modeling  of  Premixed ,  Laminar 
Steady-State  Flames .  I.  Ozone ,n  Combustion  and  Flame ,  Vol .  39,  pp.  3 01-315, 
1980. 


^ T.P .  Coffee  and  J.M.  Heimerl,  " Transport  Algorithms  for  Premixed,  Laminar 
Steady-State  Flames,”  Combustion  and  Flame,  Vol.  43,  pp.  273-289,  1981. 

.  Wamatz,  nThe  Structure  of  Laminar  Alkane-,  Alkene -,  and  Acetylene 
Flames ,n  Proceedings  of  18th  Symposium  (International)  on  Combustion , 
pp.  369-384 ,  1981. 
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II.  THE  FLAME  EQUATIONS 


In  this  section  the  equations  that  describe  a  one-dimensional,  laminar, 
premixed  flame  that  propagates  in  an  unbounded  ideal  gas  are  recalled.  The 
effects  of  radiation,  viscosity,  and  body  forces  are  ignored.  Since  the 
burning  velocity  is  small  compared  with  the  local  speed  of  sound,  the  pressure 
is  taken  to  be  constant.  The  resulting  equations4  are:* * 


Overall  Continuity: 

|£.  +  S<Pu)  , 

9 1  9x 


(1) 


Continuity  of  Species: 


9Y, 


9  Y , 


P  Tt~  +  Pu  9x~ 


9x 


(PY1V1)  +  R4M,  ,  i  -  1,  ...N  . 


i  i 


(2) 


and  Conservation  of  Energy: 
3T 


P  C  77“*  +  P  UC  r, 
p  3t  p  3x 


ST  ii*  h  ,  3(PV1V1) 
-  s?  +  i-1  hi  1 


3x 


R1M1]  . 


The  boundary  conditions  are  the  following  (t  >_  0).  For  x  =  -00; 
T  =  Tu  and  -  Yiu  (i=l,...N), 


and  for  x  =  00 ; 

ST  3Y1 


0  (i=l , • • .N) « 


(3) 


(4) 


(5) 


The  transport  properties  that  are  of  interest  are  the  diffusion 
velocities  and  the  heat  flux  q. 


III.  INPUT  COEFFICIENTS 

Thermodynamic,  kinetic  and  transport  coefficients  are  required  to  obtain 
solutions  to  the  flame  equations  of  Section  II.  The  enthalpy,  entropy  and 
heat  capacity  for  each  species  are  functions  of  temperature.  Over  the  range 
300  <_  T  <_  3000  K  they  are  adequately  described  by  the  polynominal  fits  of 
Gordon  and  Me Bride. 5 


^ F.A .  Williams ,  Combustion  Theory,  Addison-Wesley,  Reading,  MA,  1965 . 


*The  terms  are  defined  in  the  glossary . 

Gordon  and  B.J.  McBride,  " Computer  Program  for  Calculation  of  Complex 
Chemical  Equilibrium  Compositions,  Rocket  Performance,  Incident  and 
Reflected  Shocks  and  Chapman- Jouquet  Detonations, "  NASA-SP-273,  1980 
program  version . 
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The  kinetic  network,  consisting  of  14  species  and  33  reversible 
reactions,  is  taken  from  Table  1  of  Dixon- Lewis . 6  His  network  is  designed  for 
lean  to  slightly  rich  flames;  we  have  adopted  it  and  his  forward  rate 
coefficients  uncritically.  (Where  there  are  two  entries  for  forward  rate 
coefficients,  we  have  taken  those  not  enclosed  by  parentheses.)  The  rate 
coefficients  for  CHO  reacting  with  02 ,  H,  OH  and  0  are  not  listed  in  his 
table.  We  take  Dixon-Lewis ' °  suggested  values  of  3.0,  40.0,  5.0  and  10.0  x 
1012  cm3  mole-1s-1,  respectively.  Expressions  for  the  reverse  rate 
coefficients  were  computed  from  the  forward  rate  coefficients  and  the 
corresponding  equilibrium  constants.  These  constants  were  determined  from 
Gibbs  free  energy  functions. 

For  the  present  purposes  we  have  kept  the  kinetic  parameters  constant, 
i.e.,  independent  of  stoichiometry  and  of  the  nature  of  the  buffer  gas. 

To  describe  the  transport  properties  of  the  individual  species,  we  have 
adopted  the  Lennard-Jones  (or  Stockmayer)  formalism.  The  transport  parameters 
for  the  14  species  are  taken  from  Warnatz.3  The  Lennard-Jones  parameters, 

(a,  e/k),  for  argon  and  helium  are  taken  to  be  (35.42  nm,  93.3  K)  and  (25.51 
nm,  10.22  K)  ,  respectively.'7 


IV.  APPROXIMATION  METHODS 

In  this  section  we  review  the  three  approximations,  Methods  I,  V  and  VI, 
to  the  multicomponent,  polyatomic  formalism,  based  on  the  theory  of  Wang  Chang 
and  Uhlenbeck . 8-12  We  begin  with  the  most  accurate,  Method  I,  which  is  the 


^0.  Dixon-Lewis ,  "Aspects  of  the  Kinetic  Modeling  of  Methane  Oxidation  In 
Flames ,"  Proceedings  of  the  1st  Specialist  Meeting  (International)  of  the 
Combustion  Institute,  Universite '  de  Bordeaux,  France,  Paper  No.  49, 
pp.  284-289,  20-24  July  1981. 

7R.C.  Reid,  J.M.  Prausnitz,  and  T.K.  Sherwood,  The  Properties  of  Gases  and 
Liquids,  McGraw-Hill,  NY,  3rd  Ed.,  p.  678,  19?7~. 

8j.O.  Hirschf elder,  C.F.  Curtiss,  and  R.B.  Bird,  Molecular  Theory  of  Gases 
and  Liquids,  2nd  printing,  corrected  with  notes,  Jofm  WZTey  and  'Sons]  WT_ , 

tfut: - 

^ C.S .  Wang  Chang,  G.E.  Uhlenbeck,  and  J.  de  Boer,  Studies^  in  Statistical 
Mechanics,  Vol.  2,  John  Wiley  and  Sons,  NY,  1964. 

1° L.  Monchick ,  R.J.  Munn,  and  E.A.  Mason,  "Thermal  Diffusion  in  Polyatomic 
Gases:  A  Generalized  Stefan-Maxwell  Diffusion  Equation,  "  J.  Chem.  Ph]£S^_, 
Vol.  46,  pp.  3051-3058,  1966. 

H-L.  Monchick,  A.N.G.  Pereira,  and  E.A.  Mason,  "Heat  Conductivity  in 
Polyatomic  and  Polar  Gases  and  Gas  Mixtures,"  J_.  Chenr.  Phys Vol.  42, 
pp.  3241-3256,  1965. 

L.  Monchick,  K.S.  Yun,  and  E.A.  Mason,  "Formal  Kinetic  Theory  of  Transport 
Phenomena  in  Polyatomic  Gas  Mixtures,"  J.  Chem.  Phi£s .,  Vol.  39,  pp.  654-669, 
1963. 
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three-term  So nine  approximation  to  the  formalism,  Dixon-LewisA J  gives  a 
discussion  of  this  approximation,  and  we  have  generally  followed  his 
approach.  Method  V  discusses  the  case  of  constant  transport  parameters.  In 
Method  VI  we  employ  a  simplified  form  of  the  Stefan-Maxwell  equations  that 
accounts  for  molecular  diffusion  and  we  provide  a  procedure  that  accounts  for 
thermal  diffusion. 

In  each  case  we  determine  the  diffusion  velocities  V^,  i  63  1,  2,  ...  N, 
and  the  heat  flux  q. 

A.  Method  I 


1  3 

Following  Dixon-Lewis ,  J  we  can  write 


N 

q  =  X  pY 

i-1 

iVihi  - 

x  21  - 

N  RTD.T 

r  \ 

3Xi 

(6) 

o  3x 

M  Y 
i-1  i  i 

3x 

V  =  — 

1  Xi 

N  Y 

X  -1 

J-l  XJ 

D.,  I21 

ij  3x 

T 

Di  8  /, 
PY£  3x  (ln 

T). 

(7a) 

Vi  "  1  + 

Wi» 

(7b) 

a  notation  we  find  convenient.  The  diffusion  coefficients  can  be  written  in 
terms  of  the  Sonine  expansion  coefficients, 


D 


ij 


and 


D 


i 


T 


N 

l  MX 

16T  k-1  k  k  Iji 

i  25p  M.  C  iOO 
J 


M.X. 
i  i 


ai00  • 


(8) 

(9) 


Substituting  Equations  (8)  and  (9) 


into  Equation  (7)  we  find 


V . 

l 


161  £  fclji  ^  _  81  1  i_ 

25p  .  ,  v  iOO  3x  ’  5p  iOO  3x 
J  =  1 


(In  T). 


Similary,  we  can  write 

X  =  X  +  X  . 
o  o,tr  o,int 


N 
4  X 
i=l 


Xi  ailO 


-4 


N 

X 

i=l 


X.a!., 
i  lOl 


(10) 


(11) 


7  3 

G.  Dixon-Lewis,  "Flame  Structure  and  Flame  Reaction  Kinetics. 

II.  Transport  Phenomena  in  Multicomponent  Systems,"  Proc.  Roy.  Soc ., 
Vol.  A307,  pp.  111-135 ,  1968. 
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The  diffusion  velocities  and  the  heat  flux  are  now  defined  in  terms  of 
the  a^00,  a*10,  a*Q1,  and  .  If  we  use  three  terms  in  the  polynominal 

expansion,  these  are  defined  implicitly  by  the  equations 


and 


,  1 

1 

1 

1 

1 

(aioo" 

aN00’ 

a 1 1 o 

aN10’ 

a101 

=  (0  .. 

.o,x1  . 

..  xN,xL 

...  xN) 

T 

.  lhk 

lhk 

lhk 

lhk 

lhk 

(cioo  ' 

’  CN00’ 

C1 10  “* 

CN10’ 

C101 

=  (6lh 

-Lk  ‘ 

•*  6Nh  " 

6Nk’  ° 

1  \T 
aN01^ 


lhk.T 

) 


(12) 


(13) 


0)  h , k=  1 . . . N . 


The  elements  of  the  3N  by  3N  matrix  (L)  have  been  given  in  paper  I  and 
elsewhere. The  procedures  we  have  used  for  computing  these  elements  have 
been  discussed  by  us  in  paper  I. 

From  Equation  (13),  it  appears  that  we  must  solve  N  sets  of  equations  to 

obtain  the  .  in  fact,  it  has  been  recommended3’8*12  that  only  the  a 

terms  be  computed  using  this  level  of  approximation,  and  that  the  c  terms  be 
computed  using  a  one-term  expansion.  (Such  a  procedure,  in  fact,  yields  the 
Stef an-Maxwell  equations.) 


However,  this  is  not  necessary  if  we  observe  that  we  are  only  interested 
in  the  l^'s,  and  not  directly  in  the  c  terms.  We  can  convert  the  c-vector  in 
Equation  (13)  to  a  vector  whose  elements  are  in  fact  the  by  multiplying 

1  A T1  YK 

each  equation  of  (13)  by  —  and  then  summing  over  h.  After  some  algebra 

j  25p  dx 

we  find 


(L)  (1/ 


1 


••  V 


.)• 


(• 


16T  3X1 


16T  3XN 


25p  3x  ***  9  25p  3x 


,  0 


0,  0  . . .0) 

(14) 


Equation  (14)  shows  that  we  can  solve  directly  for  the  required  \J ^  and  that 
the  right-hand  side  is  not  a  function  of  either  h  or  k. 

B.  Method  V. 


This  is  the  most  simple  method  we  have  considered.  It  consists  of  an 
empirical  formula  for  the  thermal  conductivity, 

N  N 

(15) 


-ln 


A  =  0.5  £  X.A  .  +  {  E  X./A  } 

o  ..11  .,11 

i=l  i=l 

and  a  generalization  of  Fick's  Law  for  molecular  diffusion.  Specifically,  we 
have 

3  Y , 


Y.V.  =  -D.  ^ — 

li  lm  o  x 


where 


D.  = 
lm 


1  -  xi 

l  X. 

ij 


(16) 


(17) 
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In  the  Lennard-Jones  formalism  P  *  ,  where  is 

approximately  proportional  to  Since  P  «  T_b  it  is  not 

2 

unreasonable  to  assume  that  P 1. .  is  approximately  independent  of 

1 J  2 

temperature.  Generalizing,  it  is  often  assumed  that  (p  )  is  constant. 

Likewise  for  a  monatomic  gas  X  «  «p0.5/jj(2 ,2)*  ^  where  q(2,2*5*  <x  T~®  •  16  ^  14  an(j 

it  is  not  unreasonable  to  assume  that  (pX)  is  approximately  independent  of 

temperature.  We  now  outline  a  procedure  that  permits  an  a  priori  selection  of 

2 

the  constants  (pX)  and  (p  DJ  ) . 

im 


For  a  given  flame  we  know  Tu  and  Yiu,  the  temperature  and  the  mass 
fractions  of  the  unburned  mixture.  We  also  have  a  chemical  kinetics  scheme 
and  a  method  of  computing  the  specific  heats  c  ,  and  specific  enthalpies  h^. 
Since  enthalpy  is  conserved,  i.e., 


ii  Wb  -  Li  <Vi’u  ■ 


(18) 


we  can  compute  the  adiabatic  flame  temperature  Tg  by  a  numerical  trial  and 
error  procedure.  In  this  process  we  also  find  values  for  the  burned  gas  mass 
fractions,  Yig. 

At  this  level  of  approximation  we  would  like  to  use  an  appropriate 
constant  (or  "global")  c  .  So  we  will  assume  that  c  ^  =  c  =  constant,  i=l 
...  N.  Then  the  relation  h.  =  h.  +  J  _  c  .  dT  becomes  h.  =  h.  +  c 

11  To  Pi  1  i  p 

(T  -  TQ)  and  the  mixture  enthalpy  is  given  by 


N 

E 

i-1 


yihi 


N 

E 

i=l 


Yihi 


+  c  (T  -  T  ) . 
P  o 


Substituting  (19)  into  (18)  we  find 


c 

p 


llhi  (Yiu~YiB) 

T_  -  T 


(19) 


(20) 


As  a  heuristic  rule,  we  select  T  =  0.5  (Tg  +  Tu)  and  =  0.5  (Y^  +  Yiu) 
and  th^n  evaluate  and  X  using  Equations  (17)  and  (15),  respectively. 

Then  p  D,  ,  i=l  ,2  , . .  .N-l  ,  and  pX  are  evaluated.  The  diffusion  velocity  VN  is 
found  from  the  constraint 


14 A. A.  Westenbevg,  "Present  Status  of  Information  on  Transport  Properties 
Applicable  to  Combustion  Research Qnmhu&iim  and  pflam e.  Vol.  1, 
pp.  346-358 ,  1957. 
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0. 


(21) 


N 

i.i  Yivi  = 

C.  Method  VI 


For  the  thermal  conductivity  we  use  the  simplest  formula  (15)  ,  since  we 
have  found  in  paper  I  that  the  exact  choice  is  not  important.  For  the 
molecular  diffusion  velocities  l/j  we  use  the  expression  based  upon  tjie  Stefan- 
Maxwell  equations  with  the  additional  assumption  that  all  but  the  i^  species 
move  with  the  same  velocity.  The  result  is  the  formula  recommended  by 
Hirshfelder  and  Curtiss^ 


1/  .  _ 
i 


(1  -  Y±) 


3X . 
_ 1 

3  x 


(22) 


For  flames  with  low  mass  species  and  steep  temperature  gradients,  the 
neglect  of  thermal  diffusion  is  often  as  important  as  the  differences  between 
the  computational  methods.  So  we  have  generated  a  technique  that  approximates 
the  thermal  diffusion  contribution  to  the  diffusion  velocity,  W  .  We  shall 
first  derive  an  expression  for  W.  of  a  binary  mixture  and  then  generalize  the 
results. 


For  a  binary  mixture  we  identify  from  Equation  (7a) 


W 

1 


-P1  3 InT 

PYj  3x 


Since  the  thermal  diffusion  ratio  is  defined 


8 


'12 


X1X2 

PY1Y2 


T 

Di 

IT- 

12 


(23) 


(24) 


we  can  use  Equation  (22)  to  recast  Equation  (23)  in  terms  of  known  entities; 
specifically,  we  have: 

3  X 

“l  ■  k!2  %  h  »"«'  TT  ■  <25) 

Theoretical  expressions  for  have  been  derived  and  even  in  the  first 

approximation  the  expressions  are  quite  complicated.  For  the  special  case  of 
heavy  isotopes,  however,  the  first  approximation  simplifies  to° 

15(2A*  +  5)  (6C*  -  5)  (M  -  M£) 
k12  =  2A*( 1 6A*  -  12B*  +  55)  (M  +  M  ) 


X1X2 


(26) 


15J.O.  Hirsahf elder  and  C.F.  Curtiss ,  "Theory  of  Propagation  of  Flames. 

Part  I:  General  Equations,"  Proceedings  of  3rd  Symposium  (International) 
on  Combustion,  pp.  121-127,  1949. 
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Fortunately,  Equation  (26)  reproduces  thermal  diffusion  ratios  to  within 

30%,  and  so  we  use  it  as  a  simple  but  reasonable  approximation  for  k^* 
Equation  (25)  can  now  be  evaluated. 


Generalizing  Equation  (25)  we  define 


W 


i 


k 


v  I — 

im  i  3  x 


<lnT)/^r 


> 


(27) 


where ^ 


The  are  given  by  expressions  analogous  to  Equation  (26). 

From  Equation  (26)  we  can  see  that  the  influence  of  W .  will  be  the 
greater  the  larger  the  mass  differences.  Normally,  the  thermal  diffusion  ratio 
does  not  exceed  0.1  and  in  practice,  we  compute  only  for  the  lighter 
species,  i.e.,  H,  H2 ,  and  He. 

The  resulting  diffusion  velocities  ^  do  not  satisfy  Equation 

(21).  One  technique  to  satisfy  this  constraint  is  due  to  Oran  and  Boris.  ^ 
They  note  that  diffusion  velocities  enter  into  the  Stef an-Maxwell  equations  as 
differences.  They  then  conclude  that  if  a  set  of  diffusion  velocities  Vi 
satisfies  the  Stef an-Maxwell  equations,  then  so  does  the  set  (V^  +  Vc) ,  where 
Vc  is  some  constant.  The  value  of  Vc  is  chosen  such  that  the  constraint  (21) 
is  satisfied. 


ij 


(28) 


V.  RESULTS  AND  DISCUSSION 

Table  1  shows  the  computed  flame  speeds  for  mixtures  of  CH^  and  O2  in 
buffer  gases  of  N2 ,  Ar  and  He  for  the  three  transport  methods.  Quantitatively 
the  buffer  gas  is  defined  through  the  relationship 


t°2] 


[OJ  +  [X] 


=  0.21, 


(29) 


where  X  =  N2 ,  Ar  or  He  and  the  brackets  indicate  concentrations.  Only  for  X  = 
N2  in  Equation  (29)  are  stoichiometries  other  than  unity  considered.  The  span 
or  values  for  the  flame  speeds  is  seen  to  exceed  a  factor  of  seven. 


7 

°S.  Chapman  and  T.G.  Cowling,  The  t^/therngtiegl  Theorg  of  Nan^Uniformi  Gases_, 
3rd  Ed.,  Cambridge  University  Press,  1970. 

7  7 

1  E.S.  Oran  and  J.P.  Boris,  "Detailed  Modeling  of  Combustion  Systems," 
Progress  in  Energy  and  Conbustion  Sgiengg,  Vol.  7,  pp.  1-72,  1981. 
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The  largest  difference  between  Method  I,  the  most  complete  formulation  of 
the  transport,  and  the  other  two  methods  is  11%.  Compare  Method  V  and  Method 
I  for  the  case  of  stoichiometric  flame  with  N2  as  buffer.  In  paper  I,  using 
the  kinetics  given  by  Tsatsaronis , we  found  the  same  trends  for 
stoichiometric  methane-air  mixtures. 


Table  2  shows  a  comparison  between  our  results  using  the  most  precise 
transport  method  and  two  sets  of  experimental  results  that  consider  all  three 
buffer  gases.  For  stoichiometric  methane-air,  Andrews  and  Bradley,  19 
Warnatz,3  and  Dixon-Lewis  and  Islam20  report  experimental  flame  speeds  of  45, 
42.5  and  38  cm-s-1 ,  respectively. 


For  the  6.5%  and  13%  methane-air  cases  WarnatzJ  lists  experimental  values 
as  19.5  and  16 (±5)  cm-s_1.  The  reason  for  the  factor  of  two  discrepancy  in 
the  fuel-rich  case  probably  lies  in  the  fact  that  we  uncritically  adopted  a 
reaction  network  designed  for  lean  to  slightly  rich  methane/air  flames.  Other 
than  for  the  fuel-rich  case,  we  see  that  our  computed  flame  speeds  are 
realistic. 


The  fact  that  the  flame  speeds  for  the  three  transport  methods  agree  is  a 
necessary  but  not  sufficient  condition  to  judge  the  relative  effectiveness  of 
the  methods.  The  species  and  temperature  profiles  must  also  be  examined. 

Figure  1  shows  the  profiles  for  a  stoichiometric  flame  with  helium  buffer 
for  the  intermediate  species  H,  H02 ,  and  CH20.  As  can  be  seen,  the  shape  of 
each  species  profile  is  independent  of  the  transport  method.  The  differences 
in  the  peaks  are  17%,  12%,  and  11%  for  H,  H02 ,  and  CH20,  respectively.  These 
differences  are  typical  of  the  largest  found.  Most  other  species  profiles 
show  differences  less  than  these.  Indeed,  the  computed  temperature  profiles 
for  each  mixture  are  nearly  identical  for  each  of  the  three  methods. 

Figure  2  shows  that  all  three  methods  produce  nearly  identical  profiles 
for  the  intermediate  species,  H2-  Also  seen  in  Figure  2  are  the  profiles  for 
fuel  CH^ ,  and  oxidizer,  02<  Again,  aside  from  the  peaks  in  the  region  0.6  - 
1.1  mm,  the  three  transport  methods  produce  nearly  identical  profiles. 

These  peaks  were  observed  only  in  the  computations  using  helium  as  the 
buffer  gas  and  they  could  have  three  causes.  First  there  are  numerical 
errors.  These  can  be  eliminated  because,  we  routinely  increase  the  number  of 
collocation  points  by  50%,  decrease  the  time  step  error  criterion  by  a  factor 
of  three  and  increase  the  total  space  of  the  solution  interval  of  10% ,  to 
check  the  numerical  accuracy  of  the  solutions. 


18  G.  Tsatsaronis,  " Prediction  of  Propagating  Laminar  Flames  in  Methane , 
Oxygen,  Nitrogen  Mixtures,"  Combustion  and  Flame,  Vol.  33,  pp.  217-239, 
1978. 

19G.E.  Andrews  and  D.  Bradley,  "Determination  of  Burning  Velocities:  A 
Critical  Review,"  Combustion  and  Flame,  Vol.  18,  pp.  133-153,  1972. 

20G.  Dixon-Lewis  and  S.M.  Islam,  "Flame  Modeling  and  Burning  Velocity 

Measurement,"  Proceedings  of  19th  Symposium  (International)  on  Combustion, 
to  be  published. 
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Second,  there  is  a  mass-to-mole  conversion  that  can  produce  spurious 
peaks  in  the  mole  fraction  plots.  This  phenomenon  has  been  seen  and  reported 
earlier .21  It  is  briefly  discussed  in  the  appendix.  A  check  of  all  the  mass 
fraction  profiles  shows  the  same  peaks  for  O2  find  CH^  and  so  we  eliminate  this 
as  the  cause. 

Third,  there  are  thermal  diffusion  effects.  Method  I  implicitly  contains 
the  effects  of  thermal  diffusion,  Method  VI  explicitly  contains  th^m 
[Equation  (27)  ]  and  Method  V  neglects  them  entirely.  When  we  set  DRe  -  0  in 

Method  VI  the  peaks  evident  in  Figure  2  disappear  and  the  profiles  so  obtained 

for  C>2  and  CH^  resemble  those  of  Method  V.  We  conclude  that  the  peaks  in  the 
O2  and  CH^  profiles  are  due  to  thermal  diffusion  of  helium. 

As  seen  in  Section  IV,  Method  VI  contains  a  theoretically  crude 
approximation  to  thermal  diffusion.  However,  when  the  results  are  compared  to 
those  of  Method  I,  both  here  and  in  paper  I,  we  find  quite  close  agreement. 
This  indicates  that,  in  practice,  Method  VI  is  quite  precise. 

We  employ  a  relaxation  technique^  to  find  steady  state  solutions  to  the 
flame  Equations  (1-5),  and  so  we  require  trial  solutions  to  the  profiles. 
Almost  any  reasonably  well  behaved  profile  will  allow  Method  V  to  converge 
rapidly.  These  solutions  are  then  used  as  the  initial  profiles  for  Method 
VI.  In  turn,  the  solutions  for  Method  VI  are  used  as  the  trial  solutions  for 
Method  I.  Using  this  ordering  we  have  found  that  Method  VI  required  about  3 
to  5  times  more  computational  effort  than  Method  V  and  Method  I  required  10—25 
times  more.  These  values  are  for  the  studies  presented  here  and  the  ranges 
show  that  the  computational  effort  can  be  highly  problem  dependent. 

Based  upon  these  results  we  recommend  Method  V  as  the  everyday  method  of 
choice  and  Method  VI  as  the  more  precise  computational  method. 

In  order  to  emphasize  that  the  choice  of  transport  parameters  can  be 
important  we  examined  how  the  flame  speed  for  the  stoichiometric  methane-air 
case  is  affected  by  changes  in  the  transport  parameters  of  Method  V.  We 

0 

find  that  a  20%  change  in  pX  or  a  50%  change  in  p  D  will  change  the  flame 

0  02m 

speed  by  10%.  This  change  is  about  the  same  as  the  maximum  change  found 
between  methods  in  Table  1. 


J .M.  Heimerl,  " A  Contribution  to  the  Flat  Flame  Olympias :  Problem  B ,  " 
ARBRL-TR-02416,  August  1982.  See  also:  Notes  on  Numerical  Fluid  Mechanics, 
Vol.  6,  "Numerical  Methods  in  Laminar  Flame  Propagation,"  N.  Peters  and 
J.  Warmatz,  eds.,  Friedr.  Viewag  and  Son,  Braunschweig/Wiesbaden, 
pp.  71-86,  1982  (AD  A119401). 

22T.P.  Coffee  and  J.M.  Heimerl,  "Sensitivity  Analysis  For  Premixed  Laminar, 
Steady  State  Flames ,"  ARBRL-TR-02457,  January  1983  (AD  A123866 ) . 
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For  the  relative  tests  performed  here  we  have  used  the  same  sets  of 
Leonard- Jones  (or  Stockmayer)  parameters  for  the  transport,  the  same 
thermodynamic  functions  and  the  same  kinetic  network  and  parameters.  We 
varied  only  the  transport  algorithm  and  found  that  within  factors  of  10  to  20% 
Method  V  produces  accurate  results  for  both  the  profiles  and  the  flame 
speed.  Method  VI  shows  even  greater  accuracy  and  includes  the  effect  of 
thermal  diffusion.  Thus,  we  have  shown  that  the  method  used  to  approximate  the 
multicomponent,  polyatomic  transport  phenomena  for  methane-oxygen-buf f er  gas 
mixtures  is  not  critical.  In  paper  I,  we  came  to  the  same  conclusion  for  l^- 
O2-N2  mixtures.  Since  these  two  flame  types  have  widely  differing  properties 
and  since  each  flame  type  has  been  studied  over  an  extensive  range  of  gas 
mixtures,  we  infer  that  this  result  is  a  general  one. 


TABLE  1.  CALCULATED  FLAME  SPEEDS  IN  CM/S.  Stoichiometries  other 
than  unity  were  considered  only  for  the  N2  buffer. 


Buffer 

Gas 

Transport 

Method 

6.6 

%ch4 

9.5 

12  .1 

I 

17.6 

• 

00 

31.5 

V 

16.0 

39.8 

28.5 

N2 

VI 

16.7 

43.1 

30.6 

I 

78.0 

Ar 

V 

74.2 

VI 

75.0 

I 

131.7 

He 

V 

131.3 

VI 

133.0 

15 


TABLE  2.  COMPARISON  OF  OUR  COMPUTED  FLAME  SPEEDS  WITH  THOSE 
EXPERIMENTAL  VALUES  THAT  CONSIDER  ALL  THREE  BUFFER  GASES 


Ref erence/Buf f er  Gas 

N2 

Ar 

He 

Method  I 
(this  paper) 

44.8 

78.0 

131.7 

Clingman,  et  al 

40 

92 

132 

0  / 

Morgan  and  Kanez 

38.1 

72.8 

120.4 

0  7 

W.H.  Clingman,  R.S.  Brokaw,  and  R.N.  Pease ,  "Burning  Velocities  of  Methane 
with  Nitrogen-Oxygen ,  Argon-Oxygen  and  Helium-Oxygen  Mixtures ,  "  Proceedings 
of  4th  Symposium  (International)  on  Combustion,  pp.  310-313,  1953. 

Morgan  and  W.R .  Kane,  "Some  Effects  of  Inert  Diluents  on  Flame  Speeds 
and  Temperatures ,  "  Proceedings  of  4th  Symposium  (International)  on 
Combustion,  pp.  313-320,  1953. 
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x  (mm) 


Figure  1.  H,  HO2,  and  CH20  Species  Profiles  for  a  Stoichiometric  Flame 

With  Helium  Buffer  for  Transport  Method  I  (line) ,  Method  VI  (dot) , 
and  Method  V  (dash) 
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MOLE  FRACTION 


x  ( mm) 


Figure  2. 


^2  >  02’  anc*  CH4  Species  Profiles  for  a  Stoichiometric  Flame  With 
Helium  Buffer  for  Transport  Method  I  (line),  Method  VI  (dot),  and 


Method  V  (dash) 
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THE  MASS  FRACTION  TO  MOLE  FRACTION  TRANSFORMATION 
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Because  of  the  conservation  of  mass,  the  equation  describing  a  species 
profile  is  cast  in  terms  of  the  mass  fraction  for  that  species,  Y^.  The 
solutions  of  the  multispecies  equations  (and  an  energy  equation)  yield  values 
for  the  as  functions  of  space*  To  cast  the  results  into  the  more  familiar 
mole  fractions  a  transformation  +  X ^  is  used.  Herein  we  show  how  a  non¬ 
reacting  species  that  has  a  constant  mass  fraction  profile,  e.g. ,  a  diluent, 
can  have  a  mole  fraction  profile  that  is  quite  different  in  shape. 


Given  that  Y^A)  =  Y,(B)  =  where  A,  B, . . .  are  spatial  locations,  we 

ask  what  is  the  relationship  between  say,  X^(A)  and  X^(B).  In  general 


X . (A)  =  (Y  (A)/M  )/(E  Y.(A)/M.). 
i  i  i  3  J 

Then 


(A-l) 


Xi(B)  =  (Y± ( B)/M±)/(2:  Y^B)/!^) 


(A-2) 


=  (Y1(A)/Mi)/(E  YjCBj/Mj)  . 

We  require  an  expression  for  the  denominator  of  (A-2)  in  terms  of  location 
A.  Since 


and 


E  ^(Bj/M-j  =  jj*±  Yj(B)/M.  +  Yi(B)/Mi 


E  YjCAVMj  =  Y  (A)/M  +  Y±(A)/M1 


(A-3) 


(A-4) 


X1(A)/Xi(B)  =  1  +  a  , 


we  can  equate  Y^B)  and  Y^(A)  given  by  Equations  (A-3)  and  (A-4)  to  find  that 
E  Yj(B)/Mj  =  E^±  (Yj  ( B)“Yj  (A))/M  +  E  Y.(A)/M  .  (A-5) 

Substituting  (A-5)  into  (A-2)  and  using  (A-l)  we  find  that 

(A-6) 

where  a  =  E  ( (Y  .( B)-Y  . (A))/M ./(EY ,(A)/M . )) .  (A-7) 

J  Jo  J  J  J 

Note  that  if  |  a  |  «  i  then  X±(A)  -  X±  (B);  if  a  I  0  then  X±(A)  §  X±(B). 

These  last  inequalities  allow  the  possibility  that  a  flat  Y^  profile  could  be 
transformed  into  an  profile  that  is  quite  different  in  shape. 

Examples  of  this  transformation  phenomenon  have  been  found  for  N2  and  O2 
profiles  in  an  H2/02/N2  flame21  and  for  the  N20  profile  of  an  H2/N20  flame. a_1 


A-l 


G.  Dixon- Lewis,  M.M.  Sutton,  and  A.  Williams,  "Some  Reactions  of  Hydrogen 
Atoms  and  Simple  Radicals  at  High  Temperature ,"  10th  Symposium 
(International)  on  Combustion ,  pp.  495-502,  1965. 
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GLOSSARY 


•  1 

(??')*  ( 1  O* 

=  *  '  /  ft..*  J  ,  quotient  of  collision  integrals, 

ij  ij 

B*u 

=  [5ft^*  ;  -  4ft ]/^ij  *  ,  quotient  of  collision  integrals 

c*ij 

(1  2)*  (1  1)* 

=  ft  .  .  *  /ft  .  .  *  ,  quotient  of  collision  integrals. 

IJ  ij 

Cpi 

=  specific  heat,  cal-gm 

+  CP 

=  ?  ,  c  .Y.,  specific  heat  of  mixture,  cal-gm~*-K 

i=l  pi  l 

2  -1 

=  binary  diffusion  coefficient,  cm  -s 

^im 

2  -1 

=  diffusion  coefficient  of  species  i  into  a  mixture,  cm  -s 

Dij 

^iT 

2  -1 

=  multicomponent  diffusion  coefficient,  cm  -s 

=  multicomponent  thermal  diffusion  coefficient,  gm-s  ^-cm 

h± 

=  specific  enthalpy,  cal-gm 

k 

=  Boltzmann  constant  =  1.38054X10”*^  erg-K” •'■-molecule” * . 

hj 

=  thermal  diffusion  ratio. 

kim 

=  thermal  diffusion  ratio  between  species  i  and  the  rest  of  the 
mixture . 

• 

M± 

=  molecular  weight,  gm-mole 

P 

*  pressure,  atm. 

q 

—2  —  1 

=  total  energy  flux  relative  to  the  fluid  velocity,  cal-cm  -s 

Ri 

=  rate  of  production  of  ith  species  by  chemical  reaction, 
mole-cm  -s  . 

R 

=  gas  constant  =  1.9872  cal-mole  *-K  ^ . 

Ra 

3  —1—1 

=  gas  constant  =  82.05  cm  - atm- mole  -K 

t 

=  time,  s. 

s**. 

T 

-  temperature,  K. 

*»  TU 

=  temperature  of  unburned  mixture,  K. 

tb 

=  temperature  of  burned  mixture,  K. 

u 

=  fluid  velocity,  cm-s”^. 

vi 

m 

=  diffusion  velocity,  cm-s  ^ . 
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Vi 

Wi 


iu 


iB 


£±/  k 


X 

O 

X 

o,tr 

X 

o  ,int 
P 

nu.i)* 

ij 


*  diffusion  velocity  due  to  species  gradients,  cm-s"^. 

«  diffusion  velocity  due  to  the  temperature  gradient,  cm-s  * . 

»  spatial  coordinate,  cm. 

»  mole  fraction. 

33  mass  fraction. 

=  mass  fraction  of  species  i  in  the  unburned  mixture. 

«  mass  faction  of  species  i  in  the  burned  mixture. 

*  Lennard-Jones  or  Stockmayer  parameter,  K. 

*  Lennard-Jones  or  Stockmayer  parameter,  A. 

=  thermal  conductivity  of  the  mixtures,  neglecting  diffusion  effects 
cal-cm 

»  thermal  conductivity  due  to  translational  energy,  cal-cm  ^-s  *-K  * 
=  thermal  conductivity  due  to  internal  energy,  cal-cm  ^-s  ^ . 

=  fluid  density,  gm-cm  . 

(2  2 )*  (\  2)*  (l  3)* 

!  8^  •  J  ;  Qj.*  }  ;  &).*  -  collision  integral. 

5  ij  ij  ij 
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